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Abstract. The measurement of jet fragmentation functions in p+p collisions at 200 GeV is of great interest 
because it provides a baseline to study jet quenching in heavy-ion collisions. It is expected that jet quench- 
ing in nuclear matter modifies the jet energy and multiplicity distributions, as well as the jet hadrochemical 
composition. Therefore, a systematic study of the fragmentation functions for charged hadrons and iden- 
tified particles is a goal both in p-|-p and Au-I-Au collisions at RHIC. Studying fragmentation functions 
for identified particles is interesting in p-l-p by itself because it provides a test of NLO calculations at 
RHIC energies. We present a systematic comparison of jet energy spectra and fragment distributions using 
different jet-finding algorithms in p-|-p collisions in STAR. Fragmentation functions of charged and neutral 
strange particles are also reported for different jet energies. 
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1 Introduction 

RHIC experiments successfully reported the capability of 
high-pT particles to explore by parton energy loss the 
dense matter created in nuclear collisions [HE]- However, 
such measurements are affected by well-known geomet- 
] rical biases [3] and high-momentum particles may come 
from low-energy jets that did not lose much energy in the 
medium. Therefore, there has been an effort to fully re- 
construct jets in heavy-ion collisions in order to provide 
a direct measurement of the partonic kinematics, inde- 
I pendently of the fragmentation process (quenched or un- 
' quenched). The jet fragmentation function is described by 
the variable ^ = \n{Ejet/Ph)^ where Ejet is the jet energy 
and ph is the momentum of the hadron in the jet. From 
an experimental point of view we rather use, in the ab- 
sence of particle identification and at mid-rapidity, ^ = 
\n{pTjet/PT,h)^ where Prjet is the reconstructed trans- 
verse momentum of the jet. The pt spectrum of fully re- 
constructed jets in p+p collisions was measured by STAR 
with the 2003-2004 data, showing a good agreement with 
NLO calculations over the measured pr range up to 50 
GeV [HIS]. The increased statistics in the 2006 data-set 
due to the higher luminosity and detector acceptance al- 
lows one to reconstruct higher pt jets with better ac- 
curacy and to perform particle identification of jet frag- 
ments at high momentum. Current analyses [6l[7] show 
that the measurement of reconstructed jets is feasible also 
in the high-multiplicity environment of Au+Au collisions 
at RHIC. First measurements of jet fragmentation func- 
tions in STAR were done both in p+p [8l[9] and Au+Au 
[6J. In these proceedings we report a systematic study of 



fragmentation functions for charged particles in p+p for 
different jet finding algorithms. Given the inherent uncer- 
tainties in relating reconstructed jets to partons in QCD 
calculations and the complexity of jet-finding algorithms, 
it is essential to compare different algorithms to assess 
the stability of the procedures. Fragmentation functions 
of neutral strange particles {A and K^) are also reported. 
Measurements of fragmentation functions in p+p colli- 
sions are useful as a baseline for Au+Au data, in order 
to study the expected hadrochemical modification of jets 
in the medium 10 . 



2 Technical aspects 

2.1 Event selection and technical setup for Jet-Finding 

This analysis is based on the p+p data at ^^=200 GeV 
recorded by STAR during the 2005-2006 run at RHIC. 
Two online triggers are used for high pt analyses in STAR: 
the High Tower (HT) and the Jet Patch (JP) trigger. The 
HT trigger requires in addition to the minimum bias con- 
dition one tower in the Barrel Electromagnetic Calorime- 
ter (BEMC) above a transverse-energy threshold of 5.4 
GeV. The BEMC has full azimuthal (0) coverage and ex- 
tends over pseudo-rapidity — 1 < 77 < 1. Each tower has a 
size of At] X A(j) = 0.05 x 0.05. In the JP trigger the total 
energy of a larger area of Arj x Acj) = 1 x 1 is required 
to be above a threshold of 8 GeV. The total luminosities 
are 8.7 pb~^ and 11 pb~^ for JP and HT triggered events, 
respectively. 
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To avoid double-counting of electrons which leave sig- 
nals both in the Time Projection Chamber (TPC) and in 
the BEMC, we reject the BEMC energy that matches elec- 
tron candidates and keep only the corresponding track p^- 
In addition, since a charged hadron is expected to deposit 
part of its energy in the calorimeter, the energy of one 
MIP is subtracted from a tower hit by a charged hadron. 



2.2 Jet-Finding algorithms 

A jet is experimentally defined as a cluster of particles in 
the T] — (j) plane. Jet-Finding algorithms can be grouped 
into two main categories: cone and recombination algo- 
rithms (see Ref. [TT] for more details). Cone algorithms 
start from a "seed" above a given threshold, and look 
for other particles in a cone of radius R defined as = 
^J{A(jP^^rAJf^ . The seed is the energy of tracks or tow- 
ers in a given bin in the t] — (j) plane. Cone-jet finders can 
be optimized by introducing the midpoint and the split- 
ting/merging techniques [llj and iterations. When all the 
particles are used as seeds, the algorithm is called seed- 
less. An example of that is the SISCone (Seedless Infrared 
Safe Cone) algorithm p!2|[T3]. which is part of the Fast- 
Jet [H] package, a C++ code for jet reconstruction. Re- 
combination algorithms are seedless by definition and are 
in general not bound to a circular structure. The clustering 
process is based on energy /distance weighted criteria with 
a distance scale R. The "/ct" algorithm starts from merg- 
ing pairs of low-pT particles close in phase space. On the 
contrary, the "anti-/cT" starts with high momentum parti- 
cles [15 . Both "/ct" and "anti-/cT" are part of the Fast Jet 
package. The jet energy resolution can be estimated with 
a PYTHIA simulation of p+p collisions at y^=200 GeV. 
The jet finder is applied both to the PYTHIA output and 
to the output of the GEANT detector simulation plus the 
standard reconstruction chain. Fig. [T] shows the compar- 
ison between the energies of jets made of Monte Carlo 
particles (PYTHIA jets) and those built of reconstructed 
particles (RECO jets). The jet algorithm is a midpoint 
cone algorithm with R=0.7, a seed of Egeed = 0.5 GeV, a 
minimum-pT cut-off of 0.1 GeV/c and a splitting/merging 
fraction fsput = 0.5. The energy resolution is ~ 25% and 
is approximately constant for all the measured energies. 
As reported in Fig. [U the jet-energy obtained with re- 
constructed particles is on the average smaller than the 
energy of the PYTHIA-jet: this is due to the missing en- 
ergy of neutral particles like neutrons and in the re- 
constructed jets. An alternative way to estimate the jet- 
energy resolution is to study di-jets, and compare the ener- 
gies of the two reconstructed back-to-back jets. We found 
that the energy resolution calculated with this method 
is in agreement with that obtained with simulation, indi- 
cating that /cT-broadening is small compared to the jet- 
energy. The di-jet analysis is promising because it offers 
a complete model independent way to study jets. The re- 
sults shown in the following refer to jets with the highest 
energy per event. 




Fig. 1. Comparison of RECO vs PYTHIA jets: the y-axis is 
the reconstructed jet energy after GEANT and track recon- 
struction; the X-axis is the energy of jets from PYTHIA gener- 
ated particles. The vertical bars in each bin of Ejet,PYTHiA are 
the RMS's of the histograms and represent the experimental 
energy resolution. 

3 Reults 

3.1 Jet pt spectra 

We ran four different jet finders on the JP triggered events 
in order to compare their performance. The algorithms, 
described in Section \2.2\ are: (1) Leading-Order High- 
Seed Cone (LOHSCone [6l[7]) without any additional op- 
tion like splitting/merging and two choices of seeds (2 and 
4.6 GeV), (2) kr, (3) anti-/cT and (4) seedless SIS cone. 
Fig. [2] shows the jet pr spectra for the four jet finders 
with R=0.7. The jets are required to be entirely in the 
detector acceptances, therefore a cut \r]jet\ < 0.3 is ap- 
plied {r]jet is the pseudo-rapidity of the centroid of the 
jet). No correction for jet energy resolution and trigger 
bias was applied. Bin widths are chosen to be compatible 
with the experimental jet-energy resolution of ~ 25%. All 
the FastJet algorithms plus the cone with the lower seed 
{Eseed = 2 GeV) agree with each other, suggesting that 
all of them can be employed for jet analyses in p+p. The 
cone jet-finder with a high seed shows a lower yield at 
lower jet-energy because it selects jets with a high-z lead- 
ing fragment {z = PT,h/PT,jet)- This effect increases with 
the seed. 



3.2 Jet Fragmentation Functions 

3.2.1 Assessing trigger biases 

In order to assess the biases in the two online trigger selec- 
tions, the jet fragmentation functions for HT and JP data 
are compared for different jet energies (a similar study 
was done before using the neutral energy fraction of HT 
and JP jets [5j). The ^ distributions of charged particles 
are normalized to the number of jets, so that the integral 
of the spectrum is the average multiplicity of the charged 
particles in a jet of a given energy. The fragmentation 
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Fig. 2. Uncorrected reconstructed pT-spectrum of jets for dif- 
ferent jet finders: cone with Eseed = 2 GeV (blue), cone with 
Eseed = 4.6 GeV (green), kr (black), anti-/cT (red), seedless 
SIS cone (purple). 



functions are shown in Fig. [3] for JP (black circles) and 
HT (red squares) for two jet energies, 20 < Prjet < 30 
GeV and prjet > 40 GeV respectively. The bin width is 
compatible with the jet energy resolution. 
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Fig. 3. Upper plot: uncorrected ^ distribution for JP (red) 
and HT (black) triggered data for 20 < prjet < 30 GeV. The 
results are shown for the kr algorithm with R=0.7. Bottom 
plot: the same as upper plot but for prjet > 40 GeV. 



At jet energies below ~ 30 GeV the multiplicity of 
charged particles in jets is significantly different for HT 
and JP data. The observed difference is consistent with 
a stronger trigger bias in HT events: the HT trigger se- 
lects jets with a leading neutral fragment carrying a rather 
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Fig. 4. Upper plot: uncorrected ^ distributions for the four 
different jet algorithms for 20 < prjet < 30 GeV with JP 
data and R=0.7. Bottom plot: the same as upper plot but for 
PT,jet > 40 GeV. 



large fraction of the jet-energy, thus depleting the charged 
hadrons at large z. This results in a smaller multiplicity of 
charged particles for HT jets than that of JP jets (Fig. [3l 
upper panel). This effect is reduced at higher jet energies, 
since the bias is less dominant (Fig.[3l lower panel). Both 
JP and HT triggers are expected to show biases at low ^ 
(high hadron pt) which diminish for higher jet pt- JP jets 
exhibit a negligible bias for Prjet ~ 20 GeV [5j. 



3.2.2 Comparison of Jet Algorithms 



The comparison of the four jet finders at the level of frag- 
mentation functions is reported in Fig. [H for JP data and 
R=0.7. All the jet algorithms agree over the range of jet 
energies explored, from 10 GeV to above 40 GeV, therefore 
all of them can be used interchangeably for jet analyses 
in p+p. A similar study was also performed with R=0.4 
(radius used in jet reconstruction in Au+Au) with simi- 
lar performance of the jet algorithms. This suggests that 
at RHIC energies NLO effects like soft gluon radiation 
are minor, otherwise cone and recombination algorithms 
could behave differently especially in the case of smaller 
jet radii. 
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4 Fragmentation functions of identified 
particles 



STAR has excellent efficiency in finding VO particles over 
a large kinematical range {l^pr^S GeV/c), wfiicfi al- 
lows to study jet fragmentation for identified particles, 
namely A and K'^ [9j. Once tfie VO's are reconstructed 
via invariant mass analyses, tfieir decay particles are re- 
placed by tfie VO particle before tfie jet-finding step. Tfie 
jet algoritfim is a midpoint-cone witfi splitting/merging, 
= 0.5 GeV, R=0.4 and fsput = 0.5. Fig. [5] sfiows ^ 
distributions for identified A and for three jet energies. 
QCD models predict a decrease of the peak value, ^o, with 
the mass. It is found that of strange particles is smaller 
than that of charged particles (mainly pions). However, 
there are indications of an inverse mass ordering of for 
strange particles because the of A is larger than that 
of K^. Detailed systematic uncertainities associated with 
these measurements are under investigation. Similar ob- 
servations with kaons and protons were reported by other 
experiments over a large energy range p^B^ fTT j . Ongoing 
studies are focused on the comparison of fragmentation 
functions of strange particles with NLO calculations, and 
on the comparison of pr spectra of strange particles in 
jets and in the underlying event 



10<Ejet<15G&V I 



5 Summary and Outlook 



In these proceedings we report on a STAR measurement 
of uncorrected jet-pT spectra up to 50 GeV and jet frag- 
mentation functions for different jet energies in p+p at 
V^=200 GeV. Two trigger selections were compared, and 
it was shown that the High Tower and the Jet Patch trig- 
gers agree in the fragmentation functions at high jet en- 
ergies, PTjet > 40 GeV, where the trigger bias is believed 
to be gone. Despite the apparent large differences between 
different jet-finding algorithms, the measured energy spec- 
tra and fragmentation functions are very similar for differ- 
ent jet radii, suggesting that all of the algorithms can be 
interchangeably chosen in jet analyses in p+p. This is an 
important result considering p+p as baseline for jet recon- 
struction in Au+Au. In addition, the comparison of dif- 
ferent algorithms gives access to systematic errors related 
to the jet finding procedure. In the future the identifica- 
tion of jet fragments already started with strange particles 
will be pursued and extended to high-momentum parti- 
cles, exploiting the relativistic rise. The goal is to com- 
pare fragmentation functions of identified particles in jets 
in p+p and Au+Au and study the expected increase of 
hadrochemical ratios like K/ir and p/tt at high-p^ due to 
medium modifications: this would refiect in an increased 
difference between the distributions inside the jet and in 
the underlying event. Analyses of di-jets will play a signif- 
icant role in jet studies because they provide information 
on jet-energy resolution in a complete model independent 
way. 



1.2.5 



STAH PreHminary 

R<0.4 



m charged 




I 15<E_jet<20GeV 



R<0.4 



aTAR Prenminarv ^ charged 

O A (x5) 




1 O 2 23 ^ O 




Fig. 5. ^ distributions for strange particles in jets for three 
different jet energies. Gaussian fits are used to obtain the peak 
value ^0- The purple line indicates the pr-cutoff of 1 GeV for 
good Kg/ A identification. Charged hadrons are shown for com- 
parison and have a pT-cutoff of 0.5 GeV. Taken from [9 . 
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